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ABSTRACT 



6 



We present for the first time a new method to estimate the seeing using remote 
sounding from the IR night time data of the GOES 12 satellite. We discuss the derived 
correlation between the ground data and the satellite derived values from the analysis 
of the sites located at Cerro Paranal (Chile) and Roque de los Muchachos (Canary 
Islands, Spain). We get a ground-satellite correlation percentage of about 90%. Finally, 
studying the correlation between the afternoon data and the following night, we are 
able to provide a forecast of the photometric night quality. 

Key words: atmospheric effects - site testing ~ methods: statistical. 
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1 INTRODUCTION 

The capability to optimize the scientific requirements to the 
observing conditions is a challenging effort crucial to im- 
prove the performances and to increase the final efficiency of 
the system telescope-instrumentation, mainly for very large 
telescopes. The first parameter needed for this goal is the 
knowledge of the usable nights. In the last decades this 
evaluation suffered of biases due to personal judgements, 
because they were based on visual inspection. A great im- 
provement has been obtained with the use of satellite data. 
The second important parameter in the site selection and 
in the site characterization, is the image quality because, as 
well known, it affects the scientific quality of the results in 
many fields of the astronomical research. Since the first cam- 
paigns for the site selection, the criteria were based simply 
on the direct analysis of the size of the stellar images. Now, 
with the progress of the knowledge in this area, we know 
that the seeing is characterized by multiple parameters and 
affected or simply linked to several local and wide scale con- 
ditions, such as t he ext ernal air temperature and gradients 
(Lombardi et al., (|2006l ). he reafte r Paper I), pressure, wind 
velocity (Lombardi et al., (|2007l ). hereafter Paper II) and 
a link between th ese p arameters and the optical turbulence 
(Cavazzani et al. (j201l| )). It is also crucial to know the evolu- 
tion of the seeing with the time in short and long time scale, 
mainly for the future giant telescopes, for the optimization 
of the flexible scheduling. In general the testing campaigns of 
the past were expensive and time consuming and limited to 
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a few preselected sites. The use of the archive satellite data, 
instead, is of a great importance because it allows to simul- 
taneously investigate several sites on a time base of many 
years. A quantitative survey of cloud coverage and water va- 
por content above several astronomical sites have been re- 
cently obtained using both sat ellite and ground based data 
by Er asmus and van Rooyen (|2006l ). Erasmus and Sarazin 
(|2002l '). They have been among the first to demonstrate the 
capability of the satellite data t o give the amount of useful 
nights. Delia Valle et al. (|2010l . Paper III) used a similar 
analysis and, from independent data, found an agreement 
of the amount of clear nights between satellite and ground 
based data at La Palma of about 80%. An evolution of this 
analysis is presented in Cavazzani et al. ((2011, ') where we 
used a more sophisticated method and we introduced the 
concept of satellite stable night which is the best approxi- 
mation of the concept of photometric nights. In this paper 
we present for the first time an estimation of the seeing ob- 
tained using the satellite remote sounding. We analyze the 
correlation between ground based seeing and the satellite 
based seeing. This analysis is applied to two very impor- 
tant astronomical international sites such as Cerro Paranal 
(Chile) and Roque de Los Muchachos (La Palma, Canary 
Islands, Spain) in order to validate the code in two very dif- 
ferent climatic and topographic conditions. The location of 
the two sites is presented in Fig. [T] La Palma and Paranal 
are two sites in which the astronomical community built sev- 
eral facilities thanks to the good sky condition, moreover the 
community is strongly interested to maintain a high perfor- 
mance of the instrumentation. For this reason several au- 
thors focused the attention in the characterization of these 
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two sites (Murdin 1 19851 . Sarazin |2004 Varela et al. |2008| . 
ecc). 

The ESO staff was the pioneer of this topic and the long 
record of data collected at Paranal are useful tools to analyse 
the connection between astrophysical and physical environ- 
mental conditions. Differences with La Palma microclimate 
have been discussed in Paper I, Paper II, and Paper III. Pa- 
per I shows a complete analysis of the vertical temperature 
gradients and their correlation with the astronomical see- 
ing, Paper II shows an analysis of the correlation between 
wind and astronomical parameters as well as the overall long 
term weather conditions at La Palma. A statistical fraction 
of clear nights from satellite has been derived in Paper III 
using a basic approach to test the ability of the satellite to 
select clear nights. 

The mai n con ceptual difference between Erasmus & van 
Rooyen (|2006f ) analysis and Paper III is that they used 
the radio sounding vertical profile temperature as absolute 
reference to be compared with the brightness IR tempera- 
ture measured by the satellite, while we used relative devia- 
tions from the bulk of data to detect the presence of clouds. 
In particular we selected two bands sensitive to the clouds 
and plotted one band versus the other. The calibration of 
the plot gives the statistical fraction of usable nights. The 
use of the two bands separately is efficient to sense thick 
clouds, but presents some limits in case of partial coverage 
or thin clouds. For this reason we have refined the analysis 
introducing a new band sensitive to the local phenomena 
and introducing a mathematical code to correlate the three 
bands. This analysis discriminates with success changes in 
air masses showing also a first connecti on wi th seeing varia- 
tions, as presented in Cavazzani et al. (|201ll ). In this paper, 
to better analyse the correlation between satellite reflectivity 
and ground based image quality at La Palma and Paranal, 
we have used ground and satellite based data sampling the 
year 2009. We have used GOES satellite, to have homoge- 
neous results with the previous papers and easy to compare 
and discuss. Table[T]shows the geographic position and view 
angle of the satellite for each site. The paper is organized as 
follows: 

• in Section [2] we describe the used ground based data, 

• in Section [3] we describe the satellite based data, 

• in Section [4] we describe the satellite acquisition proce- 
dure, 

• in Section [5] we describe the mathematical used model, 

• in Section [6] we describe the atmospheric correlation 
function, 

• in Section [7] we describe the approach to detect small 
clouds and local perturbations, 

• in Section [H] we describe the night classification from 
satellite, 

• in Section [9] we describe the satellite seeing, 

• in Section [TO] we describe the temporal forecasting see- 
ing 

• in Section [11] we report the discussion of the results. 



2 GROUND BASED DATA 

In this analysis we have compared satellite data with the im- 
age quality in term of FWHM obtained using the differential 
image motion monitoring (DIMM) at the two sites. Data at 



Table 1. Geographic characteristics of the analyzed sites and 
GOES12 satellite. The view angle is obtained through the for- 
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Altitude 
Km 


View Angle 


Paranal 
La Palma 


-24°37' 
+28°45' 


-70°24' 
-17°52' 


2.630 
2.363 


25°00' 
64° 10' 


GOES12 


-1-0° 00' 


-75°00' 


35800 





the Observatorio del Roque de los Muchachos (ORM) are 
derived from the Robotic Differential Image Motion Mon- 
itor (known as RoboDlMM 0) of Isaac Newton Telescope 
(INT). The INT RoboDIMM, like aU classical DlMMs, re- 
lies on the method of differential image motion of telescope 
sub-apertures to calculate the seeing Fried parameter tq. 
RoboDIMM forms four separated images of the same star, 
and measures image motion in two orthogonal directions 
from which it derives four simultaneous and independent 
estimates of the seeing. The data interpretation makes use 
of the Sarazin and Roddi er's DIMM algorithm described in 
barazm & Roddier (|l990l L based on the Kolmogorov theory 
of atmospheric turbulence in the free atmosphere. At the 
present we do not have other seeing data to check possible 
local differences. But la Palma is the only site having sev- 
eral DIMM distributed along the top of the mountain. We 
are planning to follow this analysis using all the available 
DIMM data to check possible local differences for a better 
characterization of the site and to correlate direct measure- 
ments such as Cn{h) with satellite seeing. The seeing data 
at Paranal are obtained measuring the seeing of the DIMM 
at VLT observatory. The file contains also measurements of 
the flux of a reference star, in this way the flux of the star 
can trace the presence of clouds. The ground based classifi- 
cation of the night quality instead has been done using the 
night observing log of each telescope. 



3 SATELLITE BASED DATA 

In this analysis we have used GOES satellite because it is 
able to detects the infrared (IR) night time emitted radia- 
tion permitting to compare in a simultaneous way ground 
and satellite data. A detailed discussion of the p erform ance 
of this satellite is presented in Cavazzani et al. (|201ll ). The 
main advantage of GOES with respect other satellite is that 
GOES is able to observe the full Earth disk and it have 
on board an imager with five channels allowing the collec- 
tion of five simultaneous images of almost half of the Earth 
hemisphere. Our choice to use the IR channels is because 
they allow the detection of the thermal radiation emitted 
during the night from different atmospheric layers and/or 
from the soil. An appropriate choice of the wavelength al- 
lows to choose the optimal layer emission height above the 
site. If it occurs well above the soil surface, the signal be- 
comes independent of the specific soil properties and of low 
level conditions. Phenomena occurring below the selected 

^ See |http:// catserver . ing . iac. es/robodimm/ 1 
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Figure 1. Location of the two sites involved in the analysis. As seen in the inserts the selected sites presents very different topographic 
conditions; La Palma is a sharp island, Paranal is isolated peaks over the Atacama's desert. The position of GOES12 satellite projected 
on the map. Figure also shows the comparison of one matrix at Paranal and La Palma. The deformation is a result of the satellite 
observation angle. 



site (fog, low clouds...) are also avoided. In some sites, for 
example at La Palma, this aspect is of crucial importance. 



4 SATELLITE DATA ACQUISITION 

For the purposes of this work we used GOES 12 equipped 
with the imager and we have analysed the year 2009. We 
have selected the water vapour channel (B3 band) centred 
at 6.7 ij,m, able to detect high altitude cirrus clouds, the in- 
frared channel (B4 band) centred at 10.7 pm, able to detect 
middle level clouds, and the CO2 band (B6 band) centred 
at 13.3 /im, able to sense small particle such as fog, ash and 
semi-transparent high clouds. The selection of the IR chan- 
nels was done in order to detect clouds at different heights, 
because water vapor absorbs electromagnetic radiation and 
then re-emits it in various wavelength bands, in particular in 
the infrared region at 6 — 7 /im. Each obtained data is a mea- 
surements of thermal radiation emitted during the night by 



the Earth and received by the satellite detector. If clouds are 
not present, the emissions at 10.7 /im reaching the satellite 
is largely not absorbed by the atmosphere so the measured 
radiance values are due to emission from ground surface. In- 
stead when clouds are present, the emissivity drops because 
is blocked. Data are prepared by the Comprehensive Large 
Array-data Stewardship System (CLASSjj and processed 
using the free ware software McIDAS-V-1.0-beta4. For each 
site we have identified and extracted a sub-image of 1° x 1° 
having the central pixel close to the coordinates given in 
Table [T] The use of the matrix is justified by the high cor- 
relation with the single pixel. In Cavazzani et al. l|201l[ ) we 
describe this correlation, in this paper we report the corre- 
lation coefficient values for 2009: 

• Matrix correlation coefficient at Paranal= 0.97; 



2 www.class.ngdc.noaa.gov 
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Table 2. Characteristics of the G0ES12 used bands and res- 
olution at Nadir. 





Window 


Passband 

[fim] 


Resolution 
[km] 


BAND 3 


H2O 


6.50 -=-7.00 


4 


BAND 4 


IR 


10.20 ~ 11.20 


4 


BAND6 


CO2 


13.30. 


8 



• Matrix correlation coelBcient at La Palma= 0.93; 

The use of the matrix reduces the satellite noise and 
also allows us to observe the entire sky above the site. Table 
[2] shows the main characteristics of the selected bands. For 
each night we have extracted the observations at different 
hours in local time: at 17:45, 20:45, 23:45, 02:45, 05:45, 7:45, 
8:45 and 9:45. The evaluation of the amount of useful hours 
is done using all the night but 17:45 and 9:45. We have used 
both the brightness temperature at 17:45 and 9:45 to check 
a possible day-night correlation. The last column of Table [1] 
shows the satellite view angle. The insert in Figure [T] shows 
the two different projections obtained from each acquisition 
at La Palma and Paranal. 



5 THE CODE 



In Cavazzani et al. (j201ll) there is an exhaustive description 
of the mathematics approach. We summarize here the main 
definition for a complete exposition. 

The emitted monochromatic radiation intensity at a given 
A and along a vertical path at the top of the atmosphere, 
incident at a satellite instrument is given by: 



Rx = (/o)ArA {zo)+ / Ba [T{z)] Kx {z)dz 
where: 



(1) 



• K),{z) = ^^^^ ^ Weighting Function (WF) 

• B\ \T{z)] Planck function profile as function of ver- 
tical temperature profile T 

• (^o)a Emission from the earth surface at height zo 

• T\{z) Vertical transmittance from height z to space 

For a viewing path through the atmosphere at angle 9 
to the vertical, we have: 



Tx{z,e) = e 
where: 



sec6 J" K x{z)c(z) p{z)dz 



(2) 



• p{z) => Vertical Profiles of Atmospheric Density 

• K\{z) Absorption Coefficient 

• c{z) Absorbing Gas Mixing Ratio 

The emitted radiation intensities in each considered 
band A are then: 

/•oo 

-Ra3 = (/o)a3Ta3(2o) + / Bas [T{z)]Kxsiz)dz 



Rxi = (-^'o)A4-rA4(2:o) + / Ba4 [T{z)]Kx^{z)dz 



Rxe = iIo)xerxe{zo) + / Bx^ [T{z)] Kxe{z)dz 

J ZQ 

Each considered band is characterized by a weighting 
function (WF) giving the variation of the efficiency of the 
system as a function of the height. The peak of the efficiency 
specifies the layer from which the radiation is emitted and 
than the region of the atmosphere which can be sensed from 
space at fixed A. 

Assuming a standard atmosphere G0ES12 WFs assign the 
following median height values to each banc[f|: 



• BAND3: Kxs (z) 



dz 



BAND4: iCA4(^) - 
BAND6: Kx,{z) = ^-^^ 



8000m 
4000m 
3000m 



The elevation assigned by the WF depend on the lo- 
cation of the selected sites but we can assume that these 
value can be assigned to both the interested sites having an 
heights ranging between 2 and 3 Km. 



6 ATMOSPHERIC CORRELATION 
FUNCTION 

Instead to use each band separately in this analysis we have 
introduced a code to correlate the three bands. The correla- 
tion function Fc.A.(t) is given by the equation Q: 



Fc.A. = /a3 — [Ixe — Ixi] 



(3) 



In mathematical terms this model provides the bright- 
ness temperature of the B3, B4 and B6 combination, given 
by equation: 

-Ra3 + Rxa — ^^Ak 



Fc. 



- + 



t{zo) 

^ Bx,[T{z)]Kx,+ BxAT{z)]Kx, 



B,,[T{z)]K,,dz 



-(4) 



t(zo) 

The physical meaning of this model is the brightness 
temperature of the atmosphere reaching the satellite sensor 
as provided by the combination of the B3, B4 and B6 bands 
and it is given by equation Q. Therefore Fq.a., provides 
information about the atmospheric evolution of the surveyed 
site. Moreover Fc.a. provides information in both the height 
and quality of the perturbations over the surveyed site, that 
are both a function of the T brightness. In fact Figure [2] 
shows the obtained Fc.a.{^) Paranal for August 2009. Top 
Figure [5] shows the monthly trend of the three used bands. 
The central plot shows the Fc.a. of August obtained from 
the three bands, the solid gray line is the Fc.A.{t) linear 
regression. The bottom part shows the distribution of the 
clear and stable nights discussed in Section [S] 



7 DETECTION OF SUBTLE PHENOMENA 

We describe here two different algorithm introduced to de- 
tect perturbations in two cases: low level perturbations lo- 
cated spatially very close to the telescope or very far to the 



See jhttp:/ /cimss.ssec.wisc. 
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Atmospheric Stability Paranai Aug2009 




20 25 30 
Number of Events 

Figure 2. Atmosferic Correlation Function at Paranai, Au- 
gust 2009. Top figure shows the monthly plot of the three used 
bands. The central plot shows the Fc.A. of August, the solid gray 
line is the -Fc.A.(i) linear regression. The brightness tempera- 
ture is expressed in number of satellite counts as extracted with 
McIDAS-V program. The bottom part shows the distribution of 
the clear and stable nights as a function of the sensed height. 



telescope (located at the wedge of the matrix area) that 
means to have an incoming perturbation. 



7.1 Detection of small clouds in the matrix area 

In Delia Valle et al. I2OI0I the reflectivity flux has been ob- 
tained from the pixel of the matrix centred close to the co- 
ordinates of the interested site. To the aim to reduce the 
instrumental noise and to looking a wider fleld of view, we 
decided to replace the 1 pixel flux reflectivity with the mean 
value of the 1 degree matrix even centred at the coordinates 
of the interested site, moreover, to better discriminate small 
clouds distributed in the matrix area, that are missed giving 
a limitation of the model as described in the previous sec- 
tion, we computed the standard deviation of each matrix. In 
fact a high standard deviation signifies the presence of per- 
turbations in the wall area. We are able also to see incoming 
clouds approaching to the edge of the matrix area Figures [3] 
and 13] show two examples in which the average value of the 
matrix in both the figures correspond to clear nights, but the 
standard deviation of Figure [4] is high, showing a non real 
fiat distribution of the satellite counts. This is the case of 
incoming perturbation to the telescope site. Considering the 
standard deviation, we obtain the following classification: 




Figure 3. Example of a low standard deviation of pixel array. 




Figure 4. Example of a high standard deviation of pixel array. 



Standard deviation(rB) < 2a => Clear 

Standard deviation(TB) > 2a Subtle Phenomena 



7.2 Detection of local phenomena 

Finally, to better detect the presence of local phenomena 
close to the telescope, we introduce the difference between 
the mean matrix refiectivity and the single pixel refiectivity 
through the formula: 



lRs{Matrix/lPixel) = 1/^4 - hi(ipixei)\ 



(5) 



A high value of Ins shows the presence of a perturba- 
tion, even at low S/N level, not detectable using the simple 
standard deviation and the matrix average due the the high 
number of averaged pixels. 

In particular, we use this mathematics classification: 



hi 



<2a Clear 

> 2a => Subtle Phenomena 



Figure [5] represent one example of average correspond- 
ing to clear nights. In this case, the Matrix/lPixel RS 
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Figure 5. Remote Sounding between the average and the single 
pixels: example of a large difference between the average and the 
single pixels. 



shows the presence of local stationary phenomena not de- 
tected by the mean value of the matrix. By the use of both 
the Matrix /IPixel and the standard deviation of each data 
we can better to detect local phenomena and thin clouds as 
shown in Figure [6] in La Palma. In fact you see that the 
plot of B4 band evidenced by the circle is flat, typical of 
clear sky, the bottom of Figure [6] plotting Matrix /IPixel 
difference show variations indicating the presence of local 
phenomena. A check with the logbooks describe the pres- 
ence of high humidity and ice in this nights. We stress that 
these cases are rare, in fact table [3] shows the statistical re- 
sult of this analysis for the 2009 at Paranal and La Palma. 
It is given the mean monthly percentage of clear nights and 
the fraction of the clear nights with low level phenomena. 
We see that only the 1% of the 91% of clear night at Paranal 
is affected by low level phenomena, to compare with the 3% 
of the 67% at La Palma. In both cases is a very low number. 
At Paranal May shows an high number of SUBp phenom- 
ena. The check with the log gives high wind value coming 
form the see, that means high humidity justifying the high 
satellite value. Figure [9] shows the amount of clear time at 
La Palma for the 2009, in gray it is shows the percentage of 
subtle phenomena. 



8 TEMPORAL SATELLITE CLASSIFICATION 

To have reliable prediction of the time quality, we have used 
a high temporal resolution using for each night the following 
series of data: 20:45, 23:45, 2:45, 5:45, 7:45, 8:45. Using the 
brightness temperature obtained for each considered hours 
we have obtained the monthly atmospheric correlation func- 
tion. Figures [7] and [5] show the plot of the obtained temporal 
emissivity of B4 band vs B6 for the 2009 at Paranal and La 
Palma. Nights are classified according to the comment of 
the observing logs. Clear time presents high values of emis- 
sivity at both sites. As in Paper III, the classification of 
satellite time quality is done assuming that the maximum 
monthly brightness temperature in B4 band (Tg'"'^) occurs 
in clear condition. The other brightness hourly temperatures 
are correlated with Tg''^'^ when: 



Figure 6. The B4 trend (upper panel) shows no indicates while 



the RS 



Matrix /IPixel 



(bottom panel) indicates the presence of 



local phenomena; the logbooks in fact describe the presence of 
high humidity and ice. The brightness temperature is expressed in 
number of satellite counts as extracted with McIDAS-V program. 



Ground-Satellite Data Correlation 
Paranal 2009 
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Figure 7. Temporal distribution of GOES12 B4 and B6 band 
emissivity at Paranal in 2009. Sky quality classification has been 
carried out using the Paranal log. 



2a < rl^""" - Tfl < 3cr 



Clear 

Mixed 
Covered 



where Tb Brightness temperature of the 1° x 1° ma- 
trix. Table |4] shows the obtained percentage of clear, mixed 
and covered nights at Paranal and at La Palma for the 
year 2009 using all the algorithm previous described. The 
ground based classification is derived from the comments of 
the night logbook. We found a very good agreement in both 
the two sites between ground and satellite data. The last 
row of Table |4] shows the percentage of accuracy to asso- 
ciate to each obtained fraction of nights. The uncertainty is 
computed as follows: 



^Clear/ Mixed 
^Clear / Cover ed ' 
^ Mixed/ Covered 



Clear/Mixed Uncertainty 
^ Clear/Covered Uncertainty 
Mixed/Covered Uncertainty 



Usually the quality of the ground based clear nights 
is divided between photometric and spectroscopic nights. 
Also for the satellite classification we have introduced a sim- 
ilar definition introducing the concepts of stable night (pho- 



© 2009 RAS, MNRAS 000.[HfT2l 



Ground-Satellite Data Correlation 
La Palm a 2009 























J 












A** 












* 













2000 A 

4000 6000 8000 10000 12000 14000 16000 18000 



' Covered 

Clear 
> Mixed 



G0ES12 Remote Sounding 7 

Atmospheric Stability Paranal 2009 




100 150 
Number of Events 



Figure 8. Temporal distribution of G0ES12 B4 and B6 band 
emissivity at La Palma in 2009. Sky quality classification has 
been carried out using the La Palma log. 



Figure 10. Histogram of annual atmospheric stability at Paranal. 
Light-gray bars represent the stable nights, gray bars clear but 
unstable nights, black bars the nights covered. 



Table 3. Satellite Mean Monthly Percentage. Paranal and La 
Palma 2009. Subtle Phenomena (SUBp). 





Paranal 


La Palma 




Clear Time 


SUBp 


Clear Time 


SUBp 


January 


99.3 


1.0 


58.3 


3.2 


February 


99.4 


1.5 


53.7 


4.9 


March 


97.0 


1.3 


64.2 


3.2 


April 


96.8 


0.9 


82.0 


3.0 


May 


89.0 


2.1 


73.3 


2.1 


June 


80.1 


0.9 


75.9 


2.5 


July 


79.5 


1.6 


82.9 


3.3 


August 


90.5 


1.8 


86.5 


3.0 


September 


85.3 


0.3 


60.2 


2.2 


October 


85.7 


0.7 


69.2 


2.5 


November 


99.2 


0.0 


66.2 


2.8 


December 


98.3 


0.4 


32.2 


3.0 


Mean 


90.8 


1.0 


67.1 


3.1 



tometric) and clear night (spectroscopic). Considering the 



value of the Fc.A.(i) linear regression Tg' 



we define: 



Tp-n 



• \La\ < \Ip — 1 B 



< |1(t| ^ Stable 

< |2(t| ^ Clear 



• \Tb - rj''"'''""'^ > |2cr| ^ Covered 



where: 



(i) T, 



Trendline 



Brightness temperature of the Fc.a. {t) 



linear regression computed in one month 



La Palma 2009 n subtle Phenomena 3% 
■ Clear 67% 




10 11 12 



Atmospheric Stability La Palma 2009 




50 100 150 200 250 300 350 400 450 
Number of Events 

Figure 11. Histogram of annual atmospheric stability at La 
Palma. Light-gray bars represent the stable nights, gray bars clear 
but unstable nights, black bars the nights covered. 



(ii) Tp Brightness temperature of the 1° 
in one hour 



X 1° matrix 



Table [S] shows the obtained satellite mean monthly per- 
centage of clear and stable time at Paranal and la Palma. 

Figures [TOl and [TT1 show the distribution of the amount 
of clear, stable and covered time at Paranal and la Palma 
for the year 2009 obtained from the Fc.A.{t)- The ma^dmum 
of the distribution shows the sensed height and it gives the 
height in which occur the atmospheric phenomena. Figure lT^ 
shows the monthly distribution of the clear and stable nights 
at Paranal for the considered year. We see that during the 
winter months is low the percentage of stable time and the 
night is mostly only clear. 



Paranal 2009 



■ Clear 
I Stable 



Figure 9. Subtle Phenomena. La Palma 2009. 



Figure 12. Clear and stable night fractions at Paranal 2009 from 
GOES12 satellite. 
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Table 4. Temporal data analysis of Clear/Mixed/Covered time at Paranal and La Palma in 2009. 







Ground 






Satellite 






Clear 


Mixed 


Covered 


Clear 


Mixed 


Covered 


Paranal 


90.1% 


2.2% 


7.8% 


90.8% 


2.6% 


6.6% 


La Palma 


65.8% 


5.0% 


29.3% 


67.0% 


4.5% 


28.5% 






Paranal 






La Palma 




Uncertainty 


^Clear /Mixed 


^Clear / Covered 


^A'lixed/ Covered 


^Clear /Mixed 


^ Clear / C over ed 


^ Mixed/ Covered 


Percentage 


1.2% 


0.4% 


0.8% 


1.3% 


0.5% 


0.8% 



Table 5. Satellite Mean Monthly Percentage for the year 2009. 



Paranal 



La Palma 



Paranal 


Clear Time 


Stable 


Clear Time 


Stable 


January 


99.3 


83.3 


58.3 


57.6 


February 


99.4 


86.6 


53.7 


48.8 


March 


97.0 


84.4 


64.2 


57.6 


April 


96.8 


92.2 


82.0 


78.8 


May 


89.0 


79.7 


73.3 


64.0 


June 


80.1 


76.7 


75.9 


72.3 


July 


69.5 


64.9 


82.9 


75.6 


August 


90.5 


74.0 


86.5 


76.0 


September 


85.3 


70.9 


60.2 


56.8 


October 


85.7 


71.7 


69.2 


64.7 


November 


99.2 


80.8 


66.2 


59.2 


December 


98.3 


89.8 


32.2 


30.7 


Mean 


90.8 


79.6 


67.1 


61.6 


Table 6. 


Mathematical 


and statistical uncertainties of 


model in 2009 at Paranal and La Palma. 




Site 


^Totai N{G;S) 


^Statistical 




Paranal 


1.4% 


1510 


0.05% 




La Palma 


1.5% 


1510 


0.06% 





9 SATELLITE CALCULATION OF SEEING 

In Cavazzani et al. (|201lh we shown that the adopted code 
is able to discriminate with success variations of the atmo- 
spheric stability function {Fc.A.{t) ) with the optical turbu- 
lence showing the first connection between Fc.a. (t) and see- 



Table 8. Satellite FWHM at La Palma for the year 2009. We 
have calculated the correlation coefficients only for the months 
in which the RoboDIMM gives us values for more than ten 
nights. Moreover the mean seeing only refers to the clear time 
due to the fact that in the covered time the RoboDIMM does 
not work. 



Months 


bat 


Cround 


CORcoef 


January 








February 








March 


0.9 


1.0 


0.89 


April 


1.0 


1.0 


0.91 


May 


0.8 


0.8 


0.94 


June 


0.9 


0.8 


0.93 


July 


0.9 


0.9 


0.92 


August 


0.8 


0.7 


0.92 


September 


0.8 


0.7 


0.95 


October 








November 








December 









ing. In this paper we are going deeper in this analysis and, to 
better analyse the correlation between satelhte reflectivity 
and ground based image quality at La Palma and Paranal, 
we have used ground and satellite based data sampling the 
year 2009. In particular we introduce for the first time the 
concept of satellite seeing. The Fc.a. (t) measure the tem- 
perature in different atmospheric layers, and as well as the 
ground based C|. is linked to the ro and to the FWHM, it 
is possible to derive a satellite based C|., and consequently 
C^. The zero point is given empirically in this analysis. Us- 
ing the basic formulae of the seeing theory such as Fried's 
radius ro we have: 



Table 7. Satellite FWHM at Paranal for the year 2009. 



Months 


oar 


FWHM^''"-" . 

Cround 


CORcoef 


January 


0.9 


0.9 




0.91 


February 


0.8 


0.8 




0.97 


March 


0.8 


0.8 




0.88 


April 


0.7 


0.7 




0.93 


May 


0.8 


0.8 




0.93 


June 


0.8 


0.8 




0.84 


July 


0.8 


0.8 




0.79 


August 


0.9 


0.9 




0.92 


September 


0.9 


0.9 




0.96 


October 


0.8 


0.9 




0.84 


November 


0.8 


0.9 




0.95 


December 


0.8 


0.8 




0.95 



ro 



0.423 • 



1 



cos{6z 



Cr, ■ dz 



(6) 



where is the refractive index structure parameter: 



80 -lO'"^ 



The full width at half maximum is given by the follow- 
ing formula: 



FWHM = 0.1 



A 



(7) 



While the satellite FWHM is obtained through an our 
empirical model. If we assume: 

\m mTrendline\ ^2 ^2 

Us - J-B a Gt oc C„ 
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Figure 13. Ground Data-Satellite Data Correlation. Comparison 
between the FWHM calculated from the ground and the satellite 
FWHM. Paranal, January 2009 (Correlation Coefficient= 0.91). 
The satellite FWHM is calculated through the Formula |9l 



we can replace the value in the Equation[S]obtaining 
a satellite ro calculation: 

3 

I ™ ^^^^^ ^1^^ ^ I 

47r^ 



0.423 ■ ^ ■ A(6l) 
A 



I T"! rT-iTrendline I 

Ms - -is 



(8) 



Finally, using this value we get the formula for satellite 
FWHM: 



FWHMsat = 0.58-A' 



47r^ ■ A(6l) . 



(9) 



where A{6) is an empirical constant defined by the for- 



mula: 



MO)- 



IQ- 



cosO 



(10) 



where 6 is the satellite angle of view. 
Figure [TS] shows the comparison between ground based 
FWHM and satellite based FWHM computed in the same 
hours. We note the very good agreement between the two 
set of data. Figure [14] shows the dispersion of this correlation 
and its linear regression. A tentative physical interpretation 
of our correlation could be related to the Richardson number 
Ri dependent on the vertical temperature gradient. Tables 
[7] and [8] show the comparison between the seeing as given by 
the ground and those computed by satellite using the equa- 
tion [§1 We make the following observations to discuss the 
obtained values: FW H Mq^^^^^ values at Paranal are the 
DIMM data and not the VLT values; at La Palma instead 
we have calculated the correlation coefficients only for the 
months in which the RoboDIMM gives us values for more 
than ten nights. Moreover the mean seeing only refers to 
the clear time due to the fact that the RoboDIMM does not 
work during cloudy nights. 



10 TEMPORAL FORECASTING SEEING 
ANALYSIS 

In this section we analysed for the first time the possibility 
to give a forecasting value of the seeing a few hours before 
starting the observations. We have proceeded in two differ- 
ent ways to check the capability and the best procedure. 
In the first test we have correlated the brightness temper- 
ature obtained from the value at 9:45 with the brightness 
temperature obtained using the values of the nights before. 





2 
1 8 - 




16 - 
1:4 - 






I 


1:2 ■ 






Li. 


1:0 - 








08 


a 


0:6 ■ 






0:4 ■ 




0:2 ■ 




00 ■ 




0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 2,0 

Satellite FWHM 

Figure 14. Ground Data-Satellite Data Correlation. Figure 
shows the dispersion of this correlation and its linear regression. 
Paranal, January 2009 (Correlation Coefficient= 0.91). 



Paranal Jan2009 
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Days 



Figure 15. GOES 12 emissivity in B3, B4, B6 bands (upper 
panel) at Paranal for January 2009. Botton panel shows the B4, 
B6 vertical scale zoom. The brightness temperature is expressed 
in number of satellite counts as extracted with McIDAS-V pro- 
gram. 



In the second test we have correlated the brightness 
temperature obtained from the afternoon value at 17:45 and 
the same night time. Figure [15] shows the monthly distribu- 
tion of the three bands at Paranal for January 2009 (upper 
panel) and the zoom (bottom panel) in which is possible to 
see the new day point used for the forecasting seeing. We see 
that in case of hight day values the night after is stable (pho- 
tometric night). It is interesting to note that in this analysis 
we are able to give a percentage of useful nights instead of 
useful time. Table [5] shows the monthly values of the derived 
correlation at the two sites. Column 1 of tables shows the 
month, column 2 shows the number of used days, column 
3 shows the afternoon to night correlation {A ^ N is the 
correlation between the afternoon and the next night) and 
column 4 shows the night to morning correlation (A^ — >■ M is 



© 2009 RAS, MNRAS OOO.IlM 



10 S. Cavazzani et al. 



Table 9. Forecast at Paranal and at La Palma for the year 2009. A ^ N is the correlation between the afternoon 
and the next night and N ^ M \s the correlation between the morning and the night before. 



Months 


Days 


Paranal 

A N Correlation N — >■ M Correlation 


La Palma 
A — ^ A'^ Correlation N — ) 


Kl Correlation 


January 


31 


100.0 


100.0 


93.5 


96.8 


February 


28 


96.4 


100.0 


85.7 


92.9 


March 


31 


96.8 


100.0 


90.3 


93.5 


April 


30 


96.7 


100.0 


90.0 


96.7 


May 


31 


100.0 


96.8 


100.0 


96.8 


June 


30 


100.0 


93.3 


90.0 


90.0 


July 


31 


93.5 


93.5 


96.8 


96.8 


August 


31 


96.8 


96.8 


96.8 


100.0 


September 


30 


100.0 


100.0 


96.7 


100.0 


October 


31 


100.0 


100.0 


96.8 


96.8 


November 


30 


100.0 


100.0 


100.0 


100.0 


December 


31 


96.8 


96.8 


96.8 


96.8 



28DD0 



24DD0 



20DD0 



= 12DD0 



8DDD 



4D0O 



Parana! Jan2009 




Fc.a.(t) 
DIMM FWHM 



Forecast Points 



n 



Clear Time 



T 



r 



FWHM 

2.0 
1.6 
1.2 
0.8 
0.4 



1 1 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 

□ays 



Figure 16. Figure shows the trend of the Fc.A. (t)- In this function we have highlighted the new points used for the night quality forecast 
(Forecast points). Through the position of these points we can predict whether the night will be stable (Stable time) or clear (Clear 
time). In fact, Figure also shows the DIMM FWHM values at Paranal (drawn in black). We note that at these stable time points we 
have low seeing values, conversely we have high seeing values at clear time points (The ordinate on the right shows the DIMM FWHM 
values). The brightness temperature is expressed in number of satellite counts as extracted with McIDAS-V program. 



the correlation between the morning and the night before). 
For our analysis we are interested to column 3 that give the 
correlation of all the available day-night data. We see that 
at Paranal and at La Palma the correlation decrease during 
the winter months. Tables [Tol and [TTI instead show the period 
of the day in which the meteorological variation occurred. 
Column 1 shows the month, column 2 shows the percentage 
of the variation occurred in the time range between 5 p.m. 
and 6 a.m., column 3 shows the variation occurred in the 
time between 9 p.m. and 10 a.m., column 4 show the per- 
centage obtained for the not analysed day (10a.m. — 5p.m.) 
and obtained for difference. 

These numbers are obtained through the percentage of clear 
time (Table [5]) and the correlation percentages (e.g. If we 
have clear time — 70% covered time = 30%. Then we 
have A ^ N correlation — 95% and N M correlation 



= 90%. This means that 5% of the meteorological changes 
occurred between 5 p.m. and 6 a.m., 10% between 9 p.m. 
and 10 a.m. and the remaining 15% between 10 a.m. and 5 
p.m.). We see that most of the changing occur during the 
day time (from 10 a.m. to 5 p.m.) so it is possible to corre- 
late the afternoon satellite data with the next night satellite 
data. 

Figure [16] shows the trend of the Fc.A.it) obtained 
using all the brightness values and the DIMM seeing. The 
gray line is the best fit of the monthly plot. In this function 
we have highlighted the new points used for the night 
quality forecast (Forecast points). Through the position of 
these points we can predict whether the night will be stable 
(Stable time) or clear (Clear time). 

Figure [17] shows the obtained ro values from satellite. These 
values are obtained through the model described in Section 



© 2009 RAS, MNRAS OOO.nifT2l 



GOES 12 Remote Sounding 11 



Paranal Jan2009 



Table 10. Meteorological changes at Paranal for the year 



30 

25 

—20 
E 

£10 

5 



6 7 8 9 10 11 (2 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 23 29 30 31 
Days 

Figure 17. Fried's radius values calculated from satellite (ro^Sat) 
at Paranal, January 2009. The satellite FWHM calculated with 
these values has a correlation coefficient of 0.91 virith the ground 
FWHM (see Figure [T3)l . 



2009. 


Months 


5 - 6 

p.m. a.m. 


9 - 10 

p.m. a.m. 


10 - 5 

a.m. p.m. 


January 


0.0 


0.0 


100.0 


February 


100.0 


0.0 


0.0 


March 


100.0 


0.0 


0.0 


April 


100.0 


0.0 


0.0 


May 


0.0 


29.3 


70.7 


June 


0.0 


33.3 


66.7 


July 


21.5 


21.5 


57.0 


August 


35.8 


35.8 


28.3 


September 


0.0 


0.0 


100.0 


October 


0.0 


0.0 


100.0 


November 


0.0 


0.0 


100.0 


December 


50.0 


50.0 


0.0 



[9] (Equation |8]). Moreover the value of ro is computed 
taking into account the Paranal height and it is given in the 
visible range. We obtained values close to those obtained 
using ground based data. In fact, the FWHM calculated by 
the ro.Sat value has a high correlation coefficient with the 
ground FWHM (see Tables [71 E] and Figure fTil). 



11 CONCLUSION 

In this paper, as first, we have introduced the concept of 
satellite seeing using remote sounding from the IR night 
time data of the GOES 12 satellite. We have discussed the 
derived correlation between the ground data and the satel- 
lite derived values from the analysis of the sites located at 
Cerro Paranal (Chile) and Roque de Los Muchachos (Ca- 
nary Islands, Spain) for the 2009. In this analysis we used the 
Fc.A.{t) obtained correlating the monthly mean values ob- 
tained in a 1 deg matrix of each of the three selected bands. 
This functions that is a measure of the gradient of temper- 
ature among the three layers sampled by the three bands 
is, as well as the ground based C|., linked to the ro. The ro 
values derived using Fc.A.{t) (see Section [9]) at Paranal and 
La Palma are close to ground based values, in particular the 
FWHM calculated by the ro,Sat value has a high correla- 
tion coefficient with the ground FWHM (see Tables [T] [8] and 
Figure rT4l) . In this first analysis we obtained empirically the 
zero point using the DIMM seeing from each site. We have 
demonstrate that the plot of the seeing from satellite is in 
good agreement with the DIMM seeing of the same month 
(see Figure [13)) showing a correlation ranging between 80% 
and 97 % during the months at Paranal. Figure [14] shows an 
example of the dispersion of this correlation and its linear 
regression for January 2009 at Paranal. We found a better 
correlation at la Palma (89% to 95%), this is due to the fact 
that the correlation only refers to the clear time, in fact the 
RoboDIMM does not work during cloudy nights. 

Any comments we can gives about the obtained values 
due the zero point, but we have intention to refine the pro- 
cedure. As further step we are giving for the first time the 
forecasting seeing from satellite (see Figure I16p . We have 
proceeded in two ways to select the best procedure. In the 
first test we correlated the brightness temperature of the 
morning 9:45 with the values of the night before. In the 



Table 11. Meteorological changes at La Palma for the year 
2009. 



Months 


5p.m. 6a, m. 


9p.m. 10a.m. 


10a.m. 


January 


15.4 


7.7 


77.0 


February 


31.1 


15.5 


53.4 


March 


26.9 


17.9 


55.2 


April 


55.6 


18.5 


25.9 


May 


0.0 


11.9 


88.1 


June 


41.7 


41.7 


16.7 


July 


19.0 


19.0 


62.0 


August 


24.8 


0.0 


75.2 


September 


8.3 


0.0 


91.7 


October 


10.4 


10.4 


79.2 


November 


0.0 


0.0 


100.0 


December 


4.7 


4.7 


90.5 



second test we correlated the brightness temperature of the 
17:45 afternoon with the night after. The two procedures 
seems to be show similar results, with a marginal higher 
percentage in the night-morning values for both the sites, 
but for the purpose of the prediction of the image quality 
for the incoming observing night we can use the correlation 
afternoon-night. We see that at Paranal the correlation de- 
crease during the winter months, instead we found a more 
homogeneous distribution at la Palma. 

Through this afternoon-night relationship we can give an 
estimate of the photometric night quality. In fact, in Sec- 
tion [9] we have demonstrated a high correlation between the 
FWHMiiT" and the FWHM^il'Za (see Tables [3 [El and 
Figure [14]). In addition, in Section [10] we have shown how 
the afternoon data are correlated with the night data (see 
Tables]^ column 3). 

With these two results we have a model that can provide 
a satellite seeing calculation and a forecast. An interesting 
result are the values shown in Tables [TO] and 1111 The two 
tables show the monthly percentage of the changes in the 
observation conditions at the two sites during the 2009. We 
have obtained that at Paranal the variation of the meteoro- 
logical conditions occur during the day time, but the months 
of February, March and April occur in the time interval be- 
tween the 5 p.m. and 9 p.m. During the night the weather 
is almost stable. At La Palma we shown that the variation 
occur during the day. 
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